Revision of the theory of a uniform electron gas 

S. A. Mikhailov 

Institute of Physics, University of Augsburg, D-86135 Augsburg, Germany 

(Dated: April 23, 2009) 

Abstract 

The many-body theory of a uniform electron gas was developed at the end of 1950ies. The 
Coulomb interaction between electrons was considered as a perturbation and the ground state 
energy was calculated in all orders of the non-degenerate perturbation theory. We show that the 
ground state of the unperturbed system is actually strongly degenerate and that the degenerate 
perturbation theory leads to lower values of the ground state energy. 

PACS numbers: 71.10.-w, 71.10.Ca 
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The model of a uniform electron gas is very important in the many-body theory. It 
considers a system of Coulomb interacting electrons immersed in a uniformly distributed 
positive background (jellium). Although all effects related to the discrete distribution and 
the motion of atoms in the crystal lattice of a real solid are completely ignored, the electron- 
electron interaction is supposed to be treated exactly in this model. Therefore it provides 
a valuable step towards understanding of the electron-electron interaction phenomena in 
solids. The theory of the uniform electron gas was worked out by many authors more than 
50 years ago l|, y, ^, ^ and nowadays it can be found in many textbooks on the condensed 
matter theory, see e.g. Refs. 

aaa. 

In the jellium model one considers interacting electrons in the volume Q. The average 
density of the electrons = N/Q is assumed to be equal to the density of the uniform 
positive background to ensure electroneutrality. The Hamiltonian of the system reads 

2m 2 ^ |r,- — r. 



where Hq = K'y^ the kinetic energy of electrons and V = Vee+Veb+Vbb is the total Coulomb 
energy of the system which consists of the electron-electron Vee, electron-background Veb and 
background-background Vbb interaction energies. 

One of the most important physical quantity calculated by the theory is the ground state 
energy Eg of the Hamiltonian ([T]). The result of such calculations l|, the energy Eg per 
particle, measured in Rydbergs R = 

E 2 21 n Qlfi 
' - ^ + 0.0622 ln(r,) - 0.094 + . . . , (2) 



NR rl Vg 
is expressed in terms of the dimensionless density parameter 

2 / Q \ 1/3 

tq me / ^ \ ' 



(3) 



defined as the ratio of the radius of the sphere tq, which encloses one electron, to the Bohr 
radius as = fi?/me^. The r^-parameter is proportional to e^, i.e. the series ([2]) is an 
expansion in powers of the Coulomb interaction strength, and is inversely proportional to 

1/3 

Uv , i.e. the small-r^ limit corresponds to the high electron densities. The first term in the 
Eg{rs) expansion is e^-independent and corresponds to the kinetic energy of the electrons. 
The second term is the exchange Coulomb energy. All the rest terms in the expansion ([2]) 
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are referred to as the correlation energy aaa. The formula ([2]) is widely used in the 
density functional theory. 

Similar calculations were also done for a two-dimensional (2D) uniform electron gas 
In the 2D case the density parameter is defined as 

rs = ro/as = l/^/nn^aB, (4) 

where is the surface electron density, and the result reads 

Ea 1 1.2004 

^ = ^ 0.38 - 0.172r,lnr, + . . . . 5 

Equations ([2]) and ([5]) were obtained by applying the perturbation theory to the many-body 
Schrodinger equation H^! = {Hq + V)"^ = E"^, where Ho corresponded to the non-perturbed 
problem and the total Coulomb energy V was treated as a perturbation. 

In this Letter we show that the method that lead to the formulas ([2]) and was not fully 
adequate. Deriving the formulas ([2]) and ([5]) one assumed that the unperturbed ground state 
wave function is given by a single Slater determinant |\E'°) oc det |?/'k-(rj)| composed of 
the plane waves tpk.{rj) oc exp(zkj -rj) with the wavevectors k, lying under the Fermi surface 
|ki| < kp- The second (exchange) and all subsequent terms in the expansions ([2]) and ([5]) 
were obtained as the first- ((^'^|l^|\l/g)) and higher-order corrections of the non- degenerate 
perturbation theory. 

The ground state of the unperturbed Hamiltonian Hq is, however, highly degenerate. 
Consider, for example, a 2D electron gas in the shape of an L x L square and assume 
for simplicity that we have only = 4 electrons. Applying as usual the periodic boundary 
conditions we get the allowed values of the wavevector k = {2tt/ L){m, n), where m and n are 
integers. The wavevectors kj of each of these four electrons should be appropriately chosen to 
minimize the expectation value of the kinetic energy K. The first pair of electrons (with the 
spins a up and down) should evidently be placed in the state (m, n) = (0, 0). The third and 
fourth electrons have however a broad choice. The next energy "shell", with the "energy" 
rri^ + n"^ = 1, has four possible (m, ?T,)-places, {m,n) = (1,0), (0,1), (—1,0) and (0,-1), 
which give eight single-particle states (accounting for two spin projections). Two electrons 
can be distributed over the eight allowed states by = 8!/2!6! = 28 ways. All these 28 
many-body states have the same energy in the zeroth order of the perturbation theory. All 
many-body states can be classified according to their total momentum K*°* = kj, the 



total spin 5**°* and the projection of the total spin 5**°* (the conserving quantum numbers). 
Among the 28 degenerate lowest-energy states there are the states 

(a) with a finite total momentum K*°* 7^ and the finite total-spin projection S"*"* = 1 or 
-1, e.g. 

OTO OTO 000 000 
O T ' T O ' O T . O i 
000 000 OTO OiO 

[shown is the {k^, ky)-pla.ne with the unoccupied (circles) and occupied (arrows) states; 
up- and down-arrows symbolize up and down spin-polarized electrons; the symbols ltJJ. 
in the center of each drawing correspond to the always occupied state (m, n, a) = 
(0,0, T) and (0,0, i)], 

(b) with a finite total momentum K*"* 7^ and the zero spin projection 5**°* = 0, e.g. 

ono 000 000 
o M o . o M n > o M o ' 
000 000 ono 

(c) with the vanishing total momentum K*°* = and the finite total spin projection, e.g. 

OTO OiO 
O O ' T T ' O O ' 
OTO OiO 

(d) and with the zero total momentum K*°* = and the zero spin projection S**"* = 0, 
e.g. 

OTO OiO 

O O . i T > O O • 

OiO OTO 

In general, it is unknown in advance whether the ground state of the Hamiltonian H will 
have the vanishing total momentum K*°* = and the vanishing total spin 5**°* = 0, as it 
was actually assumed in ([2]), or it will correspond to a superconducting ground state 
with K*"* 7^ and/or to a partly spin-polarized ground state with 5*°* 7^ 0. The answer 
to this question can be obtained only by means of direct calculations taking into account a 
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sufficiently large number of basis many-body states (Slater determinants). Whether the true 
ground state has zero or finite K*°* and S**"* may also depend on the Coulomb interaction 



strength, i.e. on the r^-parameter, cf. Refs. 
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12|. 



But assume for the moment that the many-body ground state does correspond to K*"* = 
and 5"*°* = S*°^ = 0. The following four many-body states correspond to K*°* = and 
0: 



1^) 



QUO 



\B) 



^ o o 

i n T 



o o 

\o i o I vo T oy \o o o I 

The states {\A) + \A'))/^ and + \B'))/V2 have the total spin 5*°* 
combinations 



\B') 



^ o o 

T n i 



\o o o) 

1, while the linear 



^1 



1^) - 
V2 



1 

71 



and 



\B) - \B') 
V2 



1 

71 



o o 



o 



have the desired properties K 



tot 



LV 

0, 



i n T 



gtot ^ 



^ O i 

QUO 
\0 T 0J\ 



^ o o 



(6) 



v 



T n i 



0. The states \E'i and \l/2 have the 



(7) 



same kinetic energy, (^i|i7o|^i) = (^2|-^o|*2) = -^o = 2{h'^/2m){2n/L)^, i.e. the ground 
state of the unperturbed problem is doubly degenerate (under the conditions K*°* = 0, 
5**°* = 5**°* = 0). The solution of the Schrodinger problem should therefore be searched for 
in the form \1/ = Ci\l/i + C2'^2- The energy levels of the Hamiltonian H are then determined 
from the equation 



Eq + Vii — E V\2 
V21 Eq + V22 



E 



0, 



(8) 



where Vij = (\l/i|V^|\l/j) and Vu = V22, Vyi = V21 due to symmetry. As the off-diagonal 
Coulomb matrix elements are non-zero, in general, thus calculated ground state energy 



E„ 



Eo + Vu - \V, 
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(9) 
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is evidently lower than the one {Eq + Vii) one would get using the non-degenerate theory. 
Since both Vu and Vu are proportional to the first power of r^, the described procedure 
reduces the numerical coefficient in the exchange energy terms [0.916 in ([2]) and 1.200 in 
(^], i.e. this correction is essential even in the limit — > 0. 

In the considered simplest case of only four electrons the diagonal Coulomb matrix ele- 
ments are the same, Vn = V22, and the difference from the textbook solution ([2]), ([5]) arises 
only due to the off-diagonal Coulomb matrix elements. This is not always so. Consider for 
example another simple case with = 6 electrons. Now four electrons should be distributed 
over the eight quantum states on the last partly-occupied shell. This gives the total degen- 
eracy D of the unperturbed ground state D = qC^ = 70. Restricting ourselves by only the 
many-body states with K*°* = and 5**°* = we have to consider eight degenerate states, 
namely the states (group A) 
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The diagonal Coulomb matrix elements (n|V^|n) are the same for the states of the group A 
and for the states of the group B, but they are not equal to each other. 



Vn = V22 = V33 = Vu ^ V55 = V^66 = Vjj = Vss. (10) 

Apart from the non-zero off-diagonal terms the 8x8 Hamiltonian matrix will now have 
different diagonal terms. 

Obviously, the ground state degeneracy of the unperturbed Hamiltonian Hq dramatically 
grows with the number of particles A^, see Figure [U For example, for A^ = 410 the full 
degeneracy (including the states with the finite total momentum and spin) is D = 32C16 = 
601 080 390 and for A^ = 2042 it is D = 48C24 = 32 247603 683 100. If to take into account 
only the states with K*°* = and 5*°* = 5**°* = 0, the degeneracy will be somewhat lower 
but it will still be an exponentially huge number. 

6 



number of electrons 



FIG. 1: The ground state degeneracy of the unperturbed Hamiltonian Hq as a function of the 
particle number N. 

For some numbers (e.g. for = 10, 18, 26, 426, ...), corresponding to the fully 
occupied shells (analogs of the noble gase atoms), the non-perturbed ground state is indeed 
non-degenerate {D = 1). For these numbers the non-degenerate perturbation theory is 
formally correct. However, taking into account many-body states which include the single- 
particle states from the next empty shell, one can also reduce Eg as compared to the non- 
degenerate theory. The loss in the kinetic energy is smaller in this case than the gain in the 
Coulomb energy obtained due to the larger basis set of many-body functions. For example, 
for = 26 two-dimensional electrons (four fully occupied shells) we found that the ground 
state energy can be reduced by ~ 0.32% if to let two electrons occupy the empty fifth shell 
(these calculations have been done by the method similar to that described in ll|). 

In order to get a reasonably accurate solution of the problem at finite one had to 
diagonalize a very large Hamiltonian matrix [13]. The size of the matrix grows with the 
number of particles A^, dimensionality of space (2D 3D) and the density parameter r^. 
At ~ 1 any practically realistic calculation of this type could be made only for A^ < 10 (cf. 
Refs. jll,Q, 14, 15, 16, 17|). The impression, given by Eqs. ([2]), ([5]), that one can exactly 
solve the many-body Coulomb problem in the thermodynamic limit A^ ^ oo is evidently 
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erroneous. 

On the other hand, the uniform electron gas with a small number of particles can be 
studied in detail using the exact diagonalization technique, as has been recently done in 



few-electron parabolic quantum dots 
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11, m m 

Il7l |. One can expect, for example. 



the interaction induced Fermi liquid - Wigner solid crossover, similar to that predicted in 
quantum dots [ll, 18|. Such studies could open up new interesting directions of research 
and lead to a better understanding of the physics of Coulomb interaction in solids. 
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